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Abstract: We estimated 34 sets of Galactic model parameters for three intermediate latitude fields 
with Galactic longitudes I — 60°, / — 90°, and I = 180°, and we discussed their dependence on the 
volume. Also, we confirmed the variation of these parameters with absolute magnitude and Galactic 
longitude. The star samples in two fields are restricted with bright and unit absolute magnitude 
intervals, 4 < Mg < 5, and 5 < Mg < 6, whereas for the third field (/ — 60°) a larger absolute 
magnitude interval is adopted, 4 < Mg < 10. The limiting apparent magnitudes of star samples are 
go = 15 and go — 22.5 mag which provide space densities within distances in the line of sight ~0.9 and 
25 kpc. 

The Galactic model parameters for the thin disc are not volume dependent. However, the ones for thick 
disc and halo do show spectacular trends in their variations with volume, except for the scalelength of 
the thick disc. The local space density of the thick disc increases, whereas the scaleheight of the same 
Galactic component decreases monotonically. However, both model parameters approach asymptotic 
values at large distances. 

The relative local space density of the halo estimated by fitting the density laws to the space densities 
evaluated for all volumes is constant, except for the small ones. However it is absolute magnitude and 
Galactic longitude dependent. The axial ratio of the halo increases abruptly for the volumes where thick 
disc is dominant, whereas it approaches an asymptotic value gradually for larger volumes, indicating a 
continuous transition from disc-like structure to a spherical one at the outermost region of the Galaxy. 
The variation of the Galactic model parameters with absolute magnitude can be explained by their 
dependence on the stellar luminosity, whereas the variation with volume and Galactic longitude at 
short distances is a bias in analysis. 
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1 Introduction 



The traditional star-count analyses of the Galactic struc- 
ture have provided a picture of the basic structural and 
stellar populations of the Galaxy. E xamples and re - 
views of these analyses can be found in iBahcalllll 19861). 

iGilmore. Wvse fc Kui ikcn (1989) , 'Maiewski' ('l993 MRobin. fe j^M^frga^igl parameters. ^The results of these works 
2000|), and recently|Chcn ct al. (2001) and Siogol et all 



Our knowledge of the structure of the Galaxy, as 
deduced from star count data with colour informa- 
tion, entered now to the next step of precision with 
the advent of new surveys such as SDSS, 2MASS, DE- 
NIS, UKIDSS, VST, CFH/Megacam and Supnme. Re- 
searchers have used different methods to determine the 



The largest of the observational studies prior 
the Sloan Digital Sky Survey (SDSS) are based on pho- 
togra phic surveys. The Basle Halo Program ([Becker! 
1 19651 ) has presented the largest systematic photom etric 



survey of th e Gala x y ([d el Rio & Fcnkart 1987; Fcn kart fc Kal 
Il987l . 1 19901 . 1 19911 : IPenkartl 11989a bed). The Basle 



arc summarized in Table 1 of lKaraali. Bilir fc Hamzaoelul 
([2004). One can see that there is an improvement for 
the numerical values of the model parameters. The lo- 
cal space density and the scaleheight of the thick disc 
can be given as an example. The evaluation of the 
disc have steadily moved towards shorter scale- 
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heights, from 1.4 5 to 0.65 kpc (iGilmore fc Reidlll983l : 
IChen et al.ll200ll ) and higher local densities (2-10%). 
In many studies the range o f values for t he pa ram- 
eters is larg e . Fo r example, IChen et al.l l|200lf ) and 
Ak. Karaah fc BuseJl998l:lBuser. Rong fc Karaalil jSiegel et"alT (|200j) give 6.5-13% and 6-10%, respec- 
1999':'Karatas. Karaali fc Bus ei<200ll : lKaraah. Bilir fc^ffifel for the local space density for the thick disc. 



Halo Program photometry is currently being recal 
ibrated and reanalysed, using an improved calibra- 
tion of the RGU photometric system dBuser fc FenkartI 
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2004; Bilir, Karaah fc Busor 2004). More recent and 
future studies are being based on charge-coupled de- 
vice (CCD) survey data. 



However, one expects the most evolved numerical val- 
ues from these recent works. That is, either the range 
for this parameter should be small or a single value 
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with a small error should be given for it. It seems that 
researchers have not been able to choose the most ap- 
propriate pr ocedures for this t opic. Finally, we quote 
the work of lJuric et al.l (120051 ). Although it is a re- 
cent work, the Galactic model parameters cited in it 
are close to the ones claimed in the early works of star 
counts, i.e. 4% relative local space density and 1200 
pc scaleheight for the thick disc. However, the size of 
the field and the number of stars, 6500 deg^ and ~ 48 
million stars, are ra ther different than the ones cited 
above. Additionallv. lJuric et al.l (|2005l l admit that, fits 
applied to the entire dataset are significantly uncertain 
due to the presence of clumps and overdensities. 

Large range or different numerical values, estimated 
by different researchers, for a specific Galactic model 
parameter may be due to several reasons: 1) The Galac- 
tic model parameters are Gal actic latitude /longitude 
dependent. The two works of iBuser. Rong fc Karaalil 
(|l998l . I1999I) cited above confirm this suggestion. Al- 
though these authors give a mean value for each pa- 
rameter, there are differences between the values of a 
given para meter for different fi elds. Also, it is shown in 
the wo rks of lBilir et al.l (|2006al lbl) and lCabrera-Lavers et af] 
(j2007l ) that the Galactic model parameters are longi- 
tude dependent. 2) The Galactic model parameters 
are absolute rnagnitude (stellar luminos ity) dependent 
(jKaraah et al.l l2004l : iBilir et all l2006d ). Hence, any 
procedure which excludes this argument give Galactic 
model arguments with large ranges. 3) Distance deter- 
mination also plays an important role on the determi- 
nation of Galactic model parameters. A single colour 
magnitude diagram used for a population, for exam- 
ple, gives a mean absolute magnitude for stars with 
the same colour but with different ultraviolet excesses 
which results in a single distance for the stars in ques- 
tion. Whereas the procedure based on the colour and 
metallicity of a star gives individual absolute magni- 
tude for each star and results in more reliable distance 
determination. 

In the present study we show that additional con- 
strains are needed to be taken into consideration for 
the estimation of more reliable Galactic model param- 
eters. We will see that a model parameter is volume 
limited dependent. That is, it increases or decreases 
with the volume which covers the star sample. This 
is an explanation between the difference of numerical 
values for a specific Galactic model parameter, esti- 
mated in different works where different limiting ap- 
parent magnitudes adopted. Contrary to the expecta- 
tion and physical explanation, the relative local space 
density of the thick disc also varies with the volume. 
This is an excellent example for the bias in analysis 
which is also the main topic of this paper. 

We introduced some simplifications, i.e. we disre- 
garded the giants. However, stars in our sample are 
not brighter than go = 15 mag, hence the number of 
giants should be small, even if they exist. We quote 
iBilir et al] (|2007l ). where the percentage of giants rel- 
ative to the dwarfs are presented. Neither corrections 
for binarity and overlapping nor accounting for con- 
tamination by compact extragalactic objects could be 
made. 

In Section 2, data and reductions are presented. 



Table 1: Data for the fields. N is the number of 
stars. 
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The Galactic model parameters and their dependence 
on volume are given in Section 3. Finally, Section 4 
provides a summary and discussion. 

2 Data and reductions 

The data were taken from SDSS (DR 5) on the WE^ 
of three intermediate-latitude fields with longitudes 
I = 90° (Fl), / = 180° (F2), and I = 60° (F3). The 
Galactic latitudes of the fields Fl and F2 are equal, 
b — 50°, whereas it is 6 = 45° for F3. Also, there are 
differences between the limiting apparent magnitudes 
of the fields. For Fl and F2, 15.5 < 50 < 22.5 whereas 
for F3, 15 < 30 < 22. The field F3 has the advantage 
of providing stars with relatively small distances. Data 
for the fields are given in Table 1. The total absorption 
Am (m= u, g, r, i and z) for each band is taken from 
the query server of SDSS DR5. Thus, the de-reddened 
magnitudes, with subscript 0, are 

uo = u - Au,go = g - Ag,ro = r " Ar, (1) 
io ^i — Ai, zo ^ z — A^. 

All the colours and magnitudes mentioned here- 
after will be de-reddened. Given that the location of 
the vast majority of our targets are at distances larger 
than 0.4 kpc, it seems appropriate to apply the full 
extinction from the maps. 

According to lChen et al.l ()200lh . the distribution of 
stars in an apparent magnitude-colour diagram, go — 
(g — r)o, can be classified as follows. The blue stars 
in the range 15 < go < 18 are dominated by thick 
disc stars with turn-off at (g — r)o ~ 0.33, and for 
go > 18 the Galactic halo stars, with turn-off at (g — 
r)o ~ 0.2, become significant. Red stars, {g — r)o ~ 
1.3, are dominated by thin disc stars at all apparent 
magnitudes. 

However, the apparent magnitude-colour diagram 
and two-colour diagrams for all objects (due to short- 
age of space large amount of data not presented here) 
indicate that the stellar distributions are contaminated 
by extra-galactic objects as claimed bv lChen et al.l (j200ll 'l. 
Distinction between star/galaxy was obtained using 
command probPSF given in DR5 WEB page. There 
1 or is designated for the probability of objects be- 
ing a star or galaxy. Needless to say, separation of 
1 or strongly depends on seeing and sky bright- 
ne ss. We also applie d the "locus-projection" method 
of lJuric et al.l (2005') in order to remove hot white 
dwarfs, low-redshift quasars, and white/red dwarf un- 
resolved binaries from our sample. This procedure con- 
sists of rejecting objects at distances larger than 0.3 

^http: / /www. sdss.org/dr5/access/index. html 
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(g-r)o 



Figure 1: {g — r)o — {r — i)o two colour diagram for 
all sources. Stars lie within 0.3 mag of the locus. 



mag from the stellar locus (Figure 1). The apparent 
magnitude histogram for all objects and for stars is 
given in Figure 2. 

2.1 Absolute magnitude and distance 
determination 

We determined two sets of absolute magnitudes, 4 < 
Mg < 5 and 5 < Mg < 6, for our star sample for 
the fields Fl and F2, whereas a single but larger set 
of absolute magnitudes, 4 < Mg < 10, for F3. The 
halo and thick disc stars are dominant in the abso- 
lutely bright intervals, however the relative number 
of thin disc stars is larger in the interval 4 < Mg < 
10. We used the proc edure in lKaraali. Bilir fc Tunyell 
(|2005. hereafter KBTD where the absolute magnitude 
offset from the Hyades main sequence, AM^ , is given 
as a function of both {g — r)o colour and (So. 43 UV- 
excess, as follows: 

AMjf = C35o.43 + C2(5o.43 + Cl5o.43 + co, (2) 

where (5o.43 is the UV-excess standardized to the colour 
index {g — r)o = 0.43 in the SDSS system which corre- 
sponds to So.ao excess standardized to {B — V)o — 0.60 
in the UBV system; the coefficients Ci (i=0, 1, 2, 3) 
are functions of {g — r)o colour (Table 2) and they are 
adopted from the work of KBT, and where AM^ is 
defined as the difference in absolute magnitude of a 
program star and a Hyades star of the same {g — r)o 
colour: 

AAff = Af;-A/f . (3) 

The absolute magnitude for a Hyades star can be 
evaluated from the Hyades sequence, normalized by 
KBT (their equation 15) . This procedure which is used 




12 14 16 18 20 22 24 26 

Figure 2: Apparent magnitude histograms for all 
sources (white area) and for stars only (shaded 
area) for the fields Fl, F2, and F3. 



Table 2: Numerical values for the coefficients Ci 
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in ou r previous works (|Ak et al ] l2007al lbl: iBilir et alJ 
120071 ') has two advantages: 1) there is no need to sep- 
arate the stars into diflerent populations, and 2) the 
absolute magnitude of a star is determined by its UV- 
excess individually which provides more accurate abso- 
lute magnitudes relative to the procedure in situ where 
a specific colour magnitude diagram is used for all stars 
of the same population. When one uses the last two 
equations (2), (3) and the following one (4) which pro- 
vides absolute magnitudes for the Hyades stars it gets 
the absolute magnitude Mg of a star: 



Mg = -2.0987(5 - r) - 0.0008(ii - gf 
-h0.0842(g -r){u-g) + 7.7557(p - r) 
-0.1556(u- 3)4- 1.9714. 



(4) 



In a canonical magnitude-limited volume, the dis- 
tance to which intrinsically bright stars are visible is 
larger than the distance to which intrinsically faint 
stars are visible. The effect of this is that brighter 
stars are statistically overrepresented and the derived 
absolute magnitude s are too faint. T his effect, known 
as Malmquist bias ijMalmauistl Il92dl ) . was formalized 
into the general formula: 



Mg = Mo - (T 



dlogAjg) 
dg 



(5) 



where Mg is the assumed absolute magnitude. Mo is 
the absolute magnitude calculated for any star using 
the KBT calibration, a is the dispersion of the KBT 
calibration, and A{g) is the counts evaluated at the 
apparent magnitude go of any star. The dispersion 
in absolute magnitude calibration is about 0.25 mag 
which produces a correction of less than 0.07 mag due 
to Malmquist bias. This correction was applied to the 
SDSS data in our work. 

Combination of the absolute magnitude Mg and 
the apparent magnitude go of a star gives its distance 
r relative to the Sun, i.e. 



[g - AIg]o = 5logr - 5. (6) 

2.2 Density functions 

The density functions were evaluated for 34 volumes 
in order to estimate a set of Galactic model parame- 
ters for each sample of stars with absolute magnitudes 
4 < Mg < 5 and 5 < Mg < 6 for the fields Fl and 
F2, and with absolute magnitudes 4 < Mg < 10 for 
the field F3. The lower and upper limiting distances 
of the volumes are defined such that to obtain reli- 
able model parameters for three Galactic components, 
i.e. thin and thick discs and halo. The brighter abso- 
lute magnitude interval 4 < Mg < 5, provides larger 
distances. Hence, densities for seven volumes could 
be evaluated for each field. Whereas for the fainter 
absolute magnitude interval, 5 < Mg < 6, only six 
volumes were available. The large absolute magnitude 
interval 4 < Mg < 10 provides both short and large 
distances which results densities with eight volumes. 



There are some differences between the lower and up- 
per limiting distances of the volumes for stars with dif- 
ferent absolute magnitudes, i.e. (1.5, 3], (1.5, 5], (1.5, 
7.5], (1.5, 10], (1.5, 15], (1.5, 20], and (1.5, 25] kpc 
for 4 < Mg < 5; (1.25, 2], (1.25, 3], (1.25, 5], (1.25, 
7.5], (1.25, 10], and (1.5, 15] kpc for 5 < Mg < 6; 
and (0.9, 1.5], (0.9, 3], (0.9, 5], (0.9, 7.5], (0.9, 10], 
(0.9, 15], (0.9, 20], and (0.9, 25] kpc for 4 < Mg < 10. 
Yet, here we have not provided the corresponding ta- 
bles due to shortage of space. But all density functions 
are presented in Figure 3, as D* = \ogD + 10, where 
D = N/AVix, ^Vi,2 = {7v/180f{A/3){r^-rl); A de- 
notes the size of the field; ri and r2 denote the lower 
and upper limiting distance of the volume AVi,2; N 
is the number of stars per unit absolute magnitude; 
r* = [(rf -I- r2)/2]^/'' is the centroid distance of the 
volume AVi,2; z* — r* sin(6), b being the Galactic lat- 
itude of the field center. 



2.3 Density laws 

In this work we adopted the density laws of Basle 
Group (jBuser et al ] 119981 , 119991 ). Disc structures are 
usually parametrized in cylindrical coordinates by ra- 
dial and vertical exponentials: 

Di{R,z) = Hi exp(-|z|/^^,i) exp{~{x-Ro)/hi), (7) 

where z — zq +r sin(6), r is the distance to the object 
from the Sun, b the Galactic latitude, zq the verti- 
cal distance of the S un from the Galactic plane, 24 
pc ijJuric et al.ll2005l ). x the projection of the galac- 
tocentric distance on the Galactic plane, _ Ro th e solar 
distance from the Galactic centre (8 kpc. lReidi ri993). 
hz^i and hi are the scaleheight and scalelength, respec- 
tively, and rii is the normalized density at the solar 
radius. The suffix i takes the values 1 and 2 as long as 
the thin and thick discs are considered. 

The density law form for the spheroid component 
used in our work is as follows. 



D4R) 



exp[10.093(l - (i?/i?o)i/^)]/(i?/i?o)'/', (8) 



where R is the (uncorrected) galactocentric distance 
in spherical coordinates, and Us the normalized local 
density. R has to be corrected for the axial ratio {c/a), 



where 



R = [x^ + {z/{c/a)) 



z = r sin fo, 



2il/2 



X = [Ro + cos^ b — 2Ror cos b cos / 



1/2 



(9) 



(10) 



(11) 



with r the distance along the line of sight, and {l,b) the 
Galactic coordinates for the field under investigation. 



3 Galactic model parameters 

We estimated the local space densities and scaleheights 
for the thin and thick discs and the local space den- 
sity and the axial ratio for the halo simultaneously 



www.publish.csiro.au/journals/pasa 



5 



I 1 1 1 .1 1 r J r 1 .1 

V Fl;%Ai*;r^« 14*: 
' . J , 1 . 1 , 1 . t... ■ 


1 1 > 1 I J 1 1 1 ' 

L 'rt:1,««r£3 (-1.51 _ 

■.. 1 1 I i 1 11. !■ •■!■»■ 


' 1. ) 1 ' 1 ■ 

I>,.^ ,Pi ii^<rs2 (5,6) J 

' 1 i i 1 1 1 ' 


1 1 1 1 1 1 1 

f 1. « I -1 1 ■! '■ 



tie 4-5: 1 fiS 2 



' 1 ' 1 1 r:' (" I!' 1 ' 1 ' 
'■ 1 ■ 1 " 


' '1 ' 1 ' t ' !■< 1 ' 1 ' 

'■ 1 ■ t ' 1 ■ 1 ■ 


' I ' 1 1'' f ' I ' 
1.1,' 


' r ' T-^—r ■ 1 ' 

1 . ■ 



1 ' 1 1 -1 ' 1 1 1 ' 1 

Ft 1:5»f s» {4,q . 

!■ . 1 . i , 1 , 1 , 1 


' 1 ' 1 ' 1 ' L' 1 ' 1 
Fz ).S-tts-Wi C4B1 . 

\ 

1,1,1.1 


' 1 1 1 ' 1 ' 

F1 f;2S.«:/tSS {5v6J . 

,1.1,1. 


' 1 ' i ' 1 ' 
, ; , 1 . 1 . 



1 — ' — I — ' — I — ^ 

F1 1*«fStO (-4.6] 



: s* S s 

1 — ' — I ^ 

F2 1',S*r5*0 ,(4.5] 




a-* 

z 

Q. 

a 

6 

o. 
f 

4 

s 



F1 tiS*rs24 (4,5) Fa. i4<riiS; (4.S1 •\ Fi i^.-erslS 



_l L L 



\ F1 -lLS'«r^®5 (4,51 



,4 * MA' i* afc 5(i .6 S ,* # M, a i # -ft. » 10- .t^ 

¥1I9^' 

— 1 — — f — T — I — n 




I ' .1 ' I, ' f 

F2 1.5«-rsaS (fljij 





1 1 1 1 . 1 1 1 1 - 


' 1 ' 1 ■ 1 ' 
^''""""■~*---^^ 


' 1 ' 1 ' 1 1 .1 1 _ 


>3 osirsi Sk!>= (4,1 q"; 

1,1 


. F3 O.S<t£3kpi; (4,101 


. F3 o.a-:rs5kpe (Aio] . 


. F3 O.S ri T Skpc (4.10] . 

1,1, 



0,5 1 0-5 1 1.5 i n 1 2 3 012S4S 

I (UK) 



T T 1 J T r T [ T 1 r [ P 1 r- 

-F3 O^XfilOiipi: (4,10]- 

, , ; 1 . , . r , . , 1 , . 


■ ■ 1 1 T 1 r 1 1 i 1 


■ 1 r t [ 1 1 1 [ T T T ] 1 I 1 


■ 1 1 1 . . > 1 1 1 1 [ 1 M 1 1 1 

■ja O.S<iS2Skp6 (4,10)" 



« s * .«. * * « s <t * 12 a a i2: fei .ae 



Figure 3: Comparison of the space density functions evaluated for 34 volumes for the absolute magnitude 
intervals 4 < Mg < 5 and 5 < Mg < 6 for the fields Fl {I = 90°) and F2 (/ = 180°), and for 4 < Mg < 10 
for the field F3 (/ = 60°). 
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Table 3: The ranges of density law parameters. 
The symbol ni denotes the local space density of 
the thin disc. 



Component 


Parameter 


Unit 


Range 


Step 


Thin disc 


scalchcight 


pc 


50-450 


1 




scalclcngth 


kpc 


1-3 


0.025 


Tliick disc 


local density 


rii 


0-25 % 


0.01 




scalchcight 


pc 


350-1500 


1 




scalclcngth 


kpc 


2-5 


0.025 


Halo 


local density 


ni 


0-1 % 


0.01 




flattening 




0.1-1.0 


0.01 



by using an in situ procedure, i.e. by employing a 
method to fit the space density functions derived 
from the observations (combined for all three popu- 
lation components) with a corresponding combination 
of the adopted population-specific, analytical density 
laws. For each parameter, we determined its val- 
ues by allowing the parameter to vary within its as- 
signed range while keeping all other parameters fixed 
at their values adopted by the appropriate lowest 
model. The range of each parameter and the steps 
used in their estimation are given in Table 3. Thus, 
we have 13 sets of Galactic model parameters for thin 
and thick discs and for the halo, for each of the fields 
Fl and F2; and eight sets for F3 (Table 4). The 
sample stars with absolute magnitudes 4 < Mg < 5 
and 5 < Mg < 6 are not as close as the stars with 
4 < Mg < 10. Hence, the distance interval where the 
space densities are extrapolated to zero distance is a 
bit larger. However, they confirm the solar space den- 
sities of Htpparcos (Jahreiss & Wiclcn 1997) and they 
exhibit similar trend of stars with 4 < Mg < 10 (see 
the following sections). Hence, the model parameters 
of absolutely bright stars are as significant as the stars 
closer to the Sun. 



3.1 Galactic model parameters of the 
thin disc 

The logarithmic local space density (nl) of the thin 
disc is different for different absolute magnitudes, but 
it is constant for all sets and for three fields, for a 
given absolute magnitude interval. Numerically, n* = 
7.29 for the interval 4 < Mg < 5 and = 7.44 for 
5 < Mg < 6 and 4 < Mg < 10. However, the case is 
different for the scaleheight of the thin disc. It changes 
with the absolute magnitude and with the direction 
of the field investigated. For the absolute magnitude 
interval 4 < Mg < 5, Hi = 247 pc and Hi = 275 
pc for the fields Fl and F2, respectively. Whereas for 
5 < Mg < 6, Hi = 220 pc and Hi = 254 pc for Fl 
and F2, respectively. The scaleheight of the thin disc is 
Hi = 360 pc for the large absolute magnitude interval, 
4 < Mg < 10, for the field F3. The scalelength of the 
thin disc changes with absolute magnitude, with field 
and with volume, however no trend can be attributed 
for its variation. The most conspicuous feature for the 
scalelength of the thin disc is that it is close to 1.7 kpc 
for the large absolute magnitude interval, whereas it 
lies within 2.1 and 2.5 kpc for the brighter absolute 
magnitude intervals. 



3.2 Galactic model parameters of the 
thick disc 

The Galactic model parameters of the thick disc, i.e. 
the local space density relative to the local space den- 
sity of the thin disc (712/711), the scaleheight (i/2), and 
the scalelength (/12), vary with the location of the in- 
vestigated field, with the absolute magnitude, and with 
the volume which involves the star sample. The last 
variation could not be observed for the model param- 
eters of the thin disc. We will discuss the trend of the 
variation for each model parameter in different sec- 
tions. 

3.2.1 The local space density of the thick 
disc 

Table 5 shows that the local space density of the thick 
disc relative to the local space density of the thin disc, 
(712/711), is strongly absolute magnitude (stellar lumi- 
nosity) dependent. A slight dependence on the longi- 
tude is also conspicuous. For the field Fl (/ = 90°) 
712/711 ~ 3 and 712/11 « 10% for the absolute magni- 
tude intervals 4 < Mg < 5 and 5 < Mg < 6, respec- 
tively. For the field F2 (Z = 180°), 712/711 is rather 
close to the one for the field Fl for the brighter ab- 
solute magnitude interval, 4 < Afg < 5. However, it 
is less than the one for the fainter absolute magnitude 
interval, 712/711 ~ 7.5%. The local space density of the 
thick disc estimated for stars with 4 < Mg < 10 for 
the field F3 (Z = 60°) is nijni « 9.5%, a value which 
is generally cited in the literature. Additionally, and 
more important, is that the local space density of the 
thick disc is volume dependent. That is, it is different 
in each set in Table 5. Furthermore, there is a good 
trend in the variation of the local space density of the 
thick disc for three absolute magnitude intervals and 
for three fields. It is an increasing function of the cen- 
troid distance of the volume involving the star sample 
and its variation is steeper at short distances. This 
is more conspicuous for the field F3 where the local 
space density could be estimated at shorter distances. 
However, it approaches an asymptotic value at larger 
distances (Figure 4). 

3.2.2 The scaleheight of the thick disc 

As the scaleheight and the local space density are anti- 
correlated, the scaleheight of the thick disc (^^2) shows 
a similar trend as the inferred local space density (Fig- 
ure 5). However there are some differences between 
the trends of the scaleheights estimated for a unit ab- 
solute magnitude interval, i.e. 4 < Mg < 5 and 5 < 
Mg < 6, and for the large absolute magnitude interval, 
4 < A'lg < 10, as explained in the following. The vari- 
ation of the scaleheight for stars with 4 < Mg < 10 is 
steeper at short distances due to the reason explained 
in the Section 3.2.1. Additionally, the scaleheight for 
stars with 4 < Mg < 5 and 5 < Mg < 6 is less than 750 
pc and it approaches an asymptotic value of H2 ~ 650 
pc which is mode value of recent studies. Whereas, 
for stars with 4 < Mg < 10, the scaleheight of the 
thick disc lies between 880 and 1030 pc and one can 
not reveal any asymptotic value from this trend. 
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Table 4: 13 sets of Galactic model parameters for the fields, Field 1 and Field 2, and 8 sets for the field 
F3, as a function of absolute magnitude and volume. The columns indicate: Mg absolute magnitude 

interval, lower and upper limiting distance of the volume ri — r2 (kpc), logarithmic local space density for 
the thin disc nl, scaleheight for the thin disc Hi, scalelength for the thin disc hi, logarithmic local space 
density for the thick disc rij, scaleheight for the thick disc H2, scalelength for the thick disc /12, relative 
local space density of the thick disc n2/ni (%), logarithmic local space density for the halo rig, the axis 
ratio of the halo (c/a), relative local space density of the halo n^/ni (%), and Xmin- Scalcheights in pc. 
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The difference just claimed may be due to the slight 
difference in Galactic latitude (A6 = 5°) of F3 than the 
ones of fields Fl and F2, and shorter Galactic longitude 
{I = 60°) of the same field (see Section 4 for detail). 



3.2.3 The scalelength of the thick disc 

From Table 4, a constant value can not be attributed 
to the scalclonght (^12) of the thick disc. It varies with 
volume but the trends arc different in five panels (Fig- 
ure 6). For the sample of stars with absolute magni- 
tudes 4 < Mg < 5, it is flat, /1.2 ~ 3.75 kpc for the 
field Fl, whereas it is maximum at the least volume, 
/i2 ~ 4 kpc, and minimum at the intermediate volume, 
/i2 ~ 3 kpc for the field F2, For the fainter absolute 
magnitudes, the variation is a step function at larger 
distances and it approaches an asymptotic value. The 
scalelength is /12 = 4 kpc at two small volumes but 
/i2 = 3 kpc at four larger volumes for the field Fl, 
whereas it is almost flat, /12 ~ 4 kpc, for the field F2 
except two less values at two intermediate volumes, 
/i2 ~ 3 kpc. The scalelength of the thick disc for the 
sample of stars with 4 < Mg < 10 is less than or equal 
to the least scalelength estimated for stars in the fields 
Fl and F2, i.e. 2.3 < /12 < 3 kpc. The trend is flat 
at the intermediate distances but it decreases at short 
and large distances. 
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Figure 4: Relative local space density of the thick 
disc as a function of absolute magnitude and cen- 
troid distance for three fields. 
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Figure 5: Scaleheight of the thick disc as a function 
of absolute magnitude and centroid distance for 
three fields. 



5 
4 
3 

-o 2 

D. 

£, 5 
4 
3 
2 


5 ^ 
4 - 

'o' 

CL 

* 3 - 

2 - 

1 ^ 




F1 (4,5]" 



F2 (4,5]- 



F1 (5,6] — 



10 15 20 5 
r* (kpc) 



(d) 



F2 (5,6]- 
10 15 20 



(e) 



F3 (4,10] - 



Figure 6: Scalelength of the thick disc as a function 
of absolute magnitude and centroid distance for 
three fields. 



Table 5: Relative local space density n2/ni(%) 
and the scaleheight of the thick disc, H2 (pc) for 
the fields Fl, and F2, as a function of absolute 
magnitude, Mg, and centroid distance, r*. The 
symbol ri — r2 as in Table 4. Distances are in kpc. 
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3.3 The Galactic model parameters 
of the halo 

The Galactic model parameters of the halo, i.e. the 
local space density relative to the local space density 
of the thin disc (ns/ni) and the axial ratio (c/o), are 
also luminosity dependent and they are different for 
star samples at different locations of the Galaxy. The 
variation of (ns/ni) is the same of the variation of the 
local space density for the thick disc, i.e. it is steeper 
at short distances but, it approaches an asymptotic 
value at larger distances. This is more conspicuous for 
the field F3 where the local space density could be esti- 
mated at shorter distances (Table 4). The asymptotic 
local space densities for the fields Fl, F2 and F3 are 
0.04% (for 4 < Mg < 5), 0.15% (for 5 < Mg < 6); 
0.12% (for 4 < Mg < 5), 0.19% (for 5 < Mg < 6); 
and 0.21% (for 4 < Mg < 10) respectively. The halo is 
dominant at large distances of our Galaxy. Hence, the 
asymptotic local space densities just claimed can be at- 
tributed as the local space densities for the halo. How- 
ever, different local space densities for stars with the 
same absolute magnitudes but at different direction 
of the Galaxy, such as ns/ni = 0.04% and ns/ni = 
0.12% for the fields Fl and F2, has to be explained 
(see Section 4). 

The trend of the variation of (c/a) for five pan- 
els in Figure 7 is the same, i.e. the axial ratio of the 
halo increases monotonically up to a distance but it 
becomes fiat beyond this distance. However there are 
some differences between the variations of (c/a) for 
stars with different absolute magnitudes in different 
fields. For fields Fl and F2 where stars with a unit 
absolute magnitude interval is considered, one can ob- 
serve an asymptotic value different from each other, 
however. Whereas, the variation of (c/a) for stars in 
the field F3 gives the indication that it still increases at 
larger distances. Different (c/a) for different absolute 
magnitudes show that halo stars of different luminosity 
occupies different regions within the halo component, 
and different values for the axial ratio at different dis- 
tances show that the halo is disc-like at short distances 
but spherical at larger distances. 
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Figure 7: Axial ratio of the halo as a function of 
absolute magnitude and the centroid distance for 
three fields. 



3.4 Degeneracy 

No degeneracy could be detected in our work as tried 
to be explained in the following: 1) The variation 
of statistics for every model paxameter exhibit a 
parabola with a perfect minimum. Figure 8a gives 
the variation of for the scaleheight of thick disc 
for 4 < Mg < 5, for the field Fl, as an example. 2) 
We used the model parameters of halo estimated via 
deeper samples which get better coustrairit on the halo, 
and omitted the corresponding space density from the 
total space density of stars with 4 < Mg < 5 for the 
field Fl. We fitted the remaining space density func- 
tion to a model of two components (thin and thick 
discs) and we estimated their parameters. They are 
equal exactly the corresponding ones estimated via fit- 
ting the total density functions to a model of three 
Galactic components (Figure 8b-c). 3) Finally, we 
changed the counts of stars with 4 < Mg < 5 in Fl 
by ±\/N, adopted the parameters from the field F2 
for stars with the same absolute magnitude and we 
noticed that they do not fit (Figure 8d). The results 
obtained by the procedures applied to the mentioned 
sample of stars can be extended to other sample of 
stars. We conclude that no degeneracy exist in our 
work. 



4 Summary and Discussion 

We estimated 34 sets of Galactic model parameters 
for three fields, Fl (/ = 90°), F2 (/ = 180°), and 
F3 {I = 60°) and discussed their dependence on ab- 
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Figure 8: Four panels related to the discussion 

of degeneracy, (a) the variation of the with a 
perfect minimum, (b) comparison of the logarith- 
mic density function for stars with 4 < Mg < 5 
in the largest volume for the field Fl with best 
fitted Galactic model combined for three popula- 
tions, thin and thick discs, and halo, (c) compar- 
ison of the logarithmic density function for stars 
in panel (b) with best fitted Galactic model for 
only thin and thick discs, after omitting the space 
densities corresponding the model parameters esti- 
mated in panel (b) for the halo, (d) comparison of 
the logarithmic space density function evaluated 
by changing the star counts ±\/iV, for the stars 
with 4 < < 5 for the field Fl, with the Galactic 
model fitted to the space density function of stars 
in the field F2. Symbols (-I-) and (o) correspond 
to star counts N + y/N and N — y/N, respectively. 
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solute magnitude (stellar luminosity), Galactic longi- 
tude, and volume. The star samples in Fl and F2 are 
limited with absolute magnitudes 4 < Mg < 5 and 
5 < Mg < 6, whereas the range of the star sample in 
F3 is adopted larger, 4 < Mg < 10, to provide space 
densities at shorter distances to the Sun and to com- 
pare the resulting Galactic model parameters with the 
ones in the literature. 

All Galactic model parameters of the thin and thick 
discs, and halo are luminosity dependent. That is, 
a specific model parameter for a component of the 
Galaxy changes with the considered absolute magni- 
tude interval. The parameters are also Galactic lon- 
gitude dependent, i.e. a given parameter is differ- 
ent for the fields Fl and F2 for two samples of stars 
with the same absolute magnitudes. These findings 
are not new, only they confirm our p revious results 
(jKaraah et all |2004| : iBihr et al.l 12006^ 1 . The original 
finding in this work is the effect of the volume. Galac- 
tic model parameters varies with the volume of the 
star sample for a fixed absolute magnitude interval, 
for a field. However, the trends of the variations of the 
model parameters with absolute magnitude, location, 
and volume are different for the Galactic components, 
i.e. thin and thick discs, and halo. 



4.1 The thin disc 

For the thin disc, the logarithmic local space density 
(rii) is the same for a field for a given absolute mag- 
nitude interval, and it is constant for all volumes in 
question. Also, the scaleheight of the thin disc is con- 
stant for all volumes, however it is different for different 
fields and different absolute magnitude intervals. The 
only model parameter of the thin disc which varies 
with volume, additional to its variation with absolute 
magnitude and longitude, is the scalelength (fei). It 
lies between 2 and 2.5 kpc, without any trend how- 



4.2 The thick disc 

The behaviour of the local space density and the scale- 
height of the thick disc is different, i.e. they are limit- 
ing volume dependent, additional to absolute magni- 
tude and Galactic longitude. The local space density 
of the thick disc relative to the local space density of 
the thin disc (n2/ni) increases monotonically at short 
distances and it approaches an asymptotic value grad- 
ually. However, its range is different for different abso- 
lute magnitudes and different fields, i.e. for the fields 
Fl and F2: 2.94-3.16% and 2.74-3.22% for the abso- 
lute magnitude interval 4 < Mg < 5, 9.46-9.98% and 
6.35-7.82% for 5 < Mg < 6, respectively, and 8.63- 
9.93% for the field F3 (4 < Mg < 10). The relatively 
small nijnx values for the absolute magnitude inter- 
val 4 < Mg < 5 for the fields Fl and F2 is due to the 
reason that the thick disc stars are rary in this interval. 

The small values of the local space density remind 
us the finding of Gilmore fc Reid (1983) who claimed 
712 /ni ~ 2% for the relative l o cal sp ace density for the 
thick disc. iGilmore &: ReidI l|l983l) claimed a bright 



apparent magnitude (/ = 18 mag) and an intermedi- 
ate distance from the galactic plane in their work which 
does not contradict with the results of our work. From 
the other hand, the asymptotic value n-ijux ~ 10% for 
the fields F3 and Fl (5 < Mg < 6) is consistent with 
the work of iBensbv et all (|2005l l who estimated the 
local normalization of the thic k disc as 10% by kine- 
matical criteria. The value of ljuric et al.! (i2005| ), i.e. 
712/711 ~ 4%, estimated for a data set which covers 
6500 deg^ of the sky is close to the small values in 
our work. The distance range and the number of stars 
in their sample are 0.1-15 kpc, and ~48 million, re- 
spectively. However. lJuric et al. |(|200S) admit that fits 
applied to the entire dataset are significantly uncertain 
due to the presence of clumps and overdensities which 
is not the case in our work. 

The scaleheight of the thick disc (II2) decreases 
monotonically up to 2 ~ 4 kpc, then becomes fiat 
up to 2 ~ 15 kpc. For the fields Fl and F2, it ap- 
proaches an asymptotic value, whereas the scaleheight 
of the thick disc for the field F3 still decreases, though 
with a small gradient. The smallest and largest val- 
ues for H-i are estimated for the samples of stars with 
4 < Mg < 5 (Field F2) and stars with 4 < Mg < 10 
(Field F3), i.e. 621 and 1030 pc, respectively. As 
mentioned in the Introduction, many authors give a 
rang e for the Gal a ctic m odel parameters (see Table 
1 of iKaraali et all (|2004l )l . For example Ichen et all 
(|200ll ) give the range 580-750 pc which is rather close 
to the range of H2 in our work estimated for the unit 
absolute intervals for Fl and F2, i.e. 621-750 pc. 
From the o ther hand, the sc aleheight of the thick disc 
claimed bv iDu et al.l (I2003D . 640 pc, is equal to the 
mean of the asymptotic scaleheight estimated for the 
fields Fl and F2, i.e. 621, 624, 634, and 691 pc. 
The range of the scaleheight (880-1030 pc) estimated 



for the large absolute magnitude interval, 4 < Mg < 
10, for the field F3 is different than the ones claimed 
above. One can reveal the indication from Figure 5e 
that the asymptotic value of the H-i is lower than 880 
pc. The upper limit of the scaleheight. Hi = 1030 
pc, reminds us the values claime d for this parame- 
ter before the year 2000. Actual l y Idel Rio fc FenkartI 
(1987), Kuiikcn & Gilmord (|l989t ). lLarsenl (|l996h . and 
iBuser et al.i (|199a 119991) give Hi^l, 1, 0.98, and 0.91 
kpc, respectively. However, there are some recent works 
where l arge scaleheights w ere claimed, i.e. the upper 
limit of lSiegel et al.l (12002 ), H2 = 1 kpc, and the most 
recent work of ljuric et al (j2005l ). H2 = 1.2 kpc. 



4.3 The halo 

The behaviour of the local space density and the ax- 
ial ratio of the halo are also limiting volume, absolute 
magnitude (stellar luminosity) and Galactic longitude 
dependent. The local space density of the halo relative 
to the local space density of the thin disc (713 /ni ) is an 
increasing function for a short distance range, 2 < 2.5 
kpc for Fl and F2 and 2 < 3.5 kpc for F3, then it 
approaches an asymptotic value gradually. Halo stars 
are dominant at large distances, hence we can con- 
sider only the asymptotic local space densities. For 
the field Fl, 713/711 = 0.04% (4 < Mg < 5) is close to 
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the one claimed by iBuser etHI ()l998l . Il999l 'l. 0.0005; 
and wa/ m = 0.15% (5 < Mg £ 6) equals exa c tly to 
those o f iReid fc Maiewskil (Il993l) , iRobin et al.1 (| 19961 . 
l2000l ). and lSieeel et al.1 (|2002l ) . Additionally this yalue 
is rather close to the local space density of the halo 
claimed by many other authors. For the field F2, 
n^/ni — 0.12% (4 < Mg < 5) either equals exactly or 
is close to the local space d ensity of the halo claimed b y 
Tritton fc MortonI ll 19841). lYamaeata fc Yoshiil l| 19921 ). 
Chen et al.1 (|200ll ). and lDu et al.1 (|2003| ): and n.-^ /wi = 
0.19% (5 < Mg < 6) is yery close to corresponding lo- 
cal sp ace den sities of[G ilmore fc Reid' ('19831). iGilmorel 
(| 19841 ). and iKuiiken fc Gilmorc (1989). We should 
keep in mind that the cited local space densities for the 
halo estimated yia star count analysis refer to star sam- 
ples with a large absolute magnitude interyal, whereas 
the ones in our work are estimated for a unit absolute 
magnitude interval, i.e. 4 < Mg < 5 and 5 < Mg < 6. 
The relative local space density for the absolute mag- 
nitude interval (4 < Mg < 10) ns/m = 0.20%, is close 
to the previous one. 

The axial ratio of the halo, (c/o), show the same 
trend in five panels for three fields, however there are 
some differences between them. It increases within a 
small distance interval {z < 3.5 kpc for Fl and F2; 
and z < 5 kpc for F3) and it becomes fiat for larger 
distances. The axial ratio approaches an asymptotic 
value in panels (a)-(d), whereas there is an indication 
in panel (e) of Figure 7 that (c/a) still increases. The 
asymptotic axial ratio for the star samples with 5 < 
Mg < 6 for the fields Fl and F2, and with 4 < Mg < 5 
for the field F2 lie between 0.56 and 0.60 and it is 
close to the corresponding one cited i n recent years (c f. 
IChen et al.1 1200 ll : ISieeel et al.1 l2002l : iDu et al.1 l2003l ). 
Whereas, for stars with 4 < Mg < 5 in the field Fl, 
the asymptotic axial ratio is higher, (c/a)=0.85, and it 
is the mode value cited i n early works of star count s (cf. 
iGilmore fc ReidI (|l983l ),rKuiiken fc Gilmord (|1989| ). The 
recent works in which an axial ratio is c ited close or 
equal to the one in question are t hose of iBuser et al.l 
(|1998| . |1999| ) and lOiha et"all (| 19991 ) (their upper limit). 
Again, we should note that the cited axial ratios are 
estimated for star samples with large absolute magni- 
tude intervals, whereas our samples are restricted with 
a unit absolute magnitude interval. Hence, the agree- 
ment between our model parameters with the cited 
ones is interesting. The largest axial ratio for the star 
sample with 4 < Mg < 10 (field F3) is (c/ a) =0.73, and 
it lies in the axia l ratio interval clai med by (|Robin et al. I 
ll996l . l2000D and lOiha etaJ] (| 19991 ). 

Conclusion: The Galactic model parameters esti- 
mated for 34 sample of stars in three fields show varia- 
tion with absolute magnitude. Galactic longitude, and 
volume. The variation with absolute magnitude can 
be explained by the dependence of the model param- 
eters on stellar luminosity. Different local space den- 
sities for different absolu te magnitudes of Hipparcos 
(jJahreiss fc WielenI [l997l ) is a good confirmation for 
this argument. However, the case is different for Galac- 
tic longitude and volume. The local space density of a 
population, for example, is fixed due its definition and 
it must not change with volume or distance. Hence, 
the only explanation of its variation (in our work) can 



be done by a bias in the analysis. The gradient of 
all model parameters is large for small distances be- 
cause the contributions of three populations, i.e. thin 
and thick discs, and halo, to the space density change 
with distance. Whereas, at larger distances the con- 
tributions of (thin fc thick) discs diminish, and the 
gradient approaches zero. Thus the asymptotic value 
of the corresponding variation can be adopted as the 
Galactic model parameter in question. 

The perfect space density functions in our work re- 
ject any possible presence of clumps and overdensities 
in the fields investigated. Also, the variation of ior 
a given model parameter can be fitted to a parabola 
which provides a perfect minimum, i.e. Xmm- This is 
a clue for the absence of any degeneracy in our work. 
Also, the combination of two fields confirm the non- 
degeneracy. 
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